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Abstract 

Tfie effect of an externally applied flow on symmetry degenerated waves propagating into opposite 
directions and standing waves that exchange stability with the traveling waves via mixed states is 
analyzed. Wave structures that consist of spiral vortices in the counter rotating Taylor-Couette 
system are investigated by full numerical simulations and explained quantitatively by amplitude 
equations containing quintic coupling terms. The latter are appropriate to describe the influence of 
inversion symmetry breaking perturbations on many oscillatory instabilities with 0(2) symmetry. 
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Many nonlinear structure forming systems that are driven out of equilibrium show a 
transition to traveling waves (TWs) as a result of an oscillatory instability ij]. In the 
presence of inversion symmetry in one or two spatial directions also a standing wave (SW) 
solution bifurcates in addition to the symmetry degenerate, oppositely propagating TWs 
at the same threshold. Moreover, depending on the parameters one can have a stability 
exchange between TWs and SWs as the driving rate varies. The stability transfer is mediated 
by mixed patterns that establish in the solution space a connection between a pure TW and 
a pure SW. 

Here we investigate how externally tunable symmetry breaking perturbations change 
stability, bifurcation properties, and the spatiotemporal behavior of the afore mentioned 
structures. To be concrete we investigate vortex waves in the annular gap between counter 
rotating concentric cylinders of the Taylor-Couette system . To that end we performed 
full numerical simulations of the Navier-Stokes equations with methods described in [4]. We 
elucidate that and how these results can be explained quantitatively by amplitude equations 
which contain only quintic order terms. 

The perturbation is realized in our system by an externally enforced axial through-flow 
that can easily be controlled experimentally. However, our results concerning the influence 
of inversion symmetry breaking perturbations on TWs, SWs, and on the mixed states are 
more general: the quintic order amplitude equations with small symmetry breaking terms 
apply to all kinds of 0(2) symmetric oscillatory instabilities in the presence of inversion 
symmetry breaking perturbations. 

Structures - Without symmetry breaking through-flow the following oscillatory vortex 
structures occur at small driving Forward bifurcating TWs consisting of left 

handed spiral vortices (L-SPI) or of right handed spiral vortices (R-SPI) that are mirror 
images of each other. L-SPI (R-SPI) travel in the annulus between the two cylinders into 
(opposite to) the direction of the rotation frequency vector of the inner one, i.e., in our 
notation upwards (downwards) [J], (ii) Forward bifurcating SWs that consist of an equal- 
amplitude nonlinear combination of L-SPI and R-SPI that are called ribbons (RIBs) in the 

nn 

Taylor-Couette literature [7|, |8|. (iii) So-called cross-spirals (CR-SPI), i.e., combinations of 
L-SPI and R-SPI with different amplitudes. They provide a stability transferring connection 
between TW and SW solution branches 6]. The vortex structures are axially and 

azimuthally periodic with wave numbers k = 2tt/X and M, respectively, with A = 1.3 and 
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M = 2 throughout this paper. They rotate with characteristic constant angular velocities 
as a whole into the same direction as the inner cylinder 1^. Thereby they are forced to 
propagate axially with the exception of RIB vortices which rotate only but do not propagate. 

Order parameters - Figures show bifurcation diagrams of SPI, CR-SPI, and RIB 
solutions in a system of radius ratio rj = 1/2 for fixed Ri = 240 versus the reduced distance 
/i = (-R2 ~ R2) / 1-^21 from the spiral onset i?2 absence of through-flow, Re = 0. Here, 

Ri and R2 are the Reynolds numbers defined by the rotational velocities of the inner and 
outer cylinder, respectively, and Re is the Reynolds number of the imposed axial through- 
flow. In Figs. [H and Ej the influence of a small through-flow {Re = 0.02) is compared with 
the situation without through-flow. Order parameters in Fig. [T]are the squared amplitudes 
of the dominant critical modes 

«2,i(t) = lAle-^"^-**, U2,-i{t) = iBle-''^^' (1) 

in the double Fourier decomposition of the radial velocity u at mid-gap in azimuthal and 
axial direction. The indices 2 and ±1 identify azimuthal and axial modes, respectively. The 
linear stability analysis of the basic state consisting of a superposition of circular Couette 
flow and of annular Poiseuille flow in axial direction shows that the growth rates of these 
modes become positive at the respective bifurcation thresholds fiA and /i^ of M = 2 L-SPI 
and R-SPI, respectively. Note also that for all relaxed vortex structures investigated here 
the moduli and frequencies of U2,±i in Eq. ([T]) are constant. 

In Fig. [2] we show in addition the bifurcation diagrams of the combinations 

S = {\A\^ + \B\^) /2, D = {\A\^ - \B\^) /2 (2) 

since they are convenient to describe in particular CR-SPI. 

Bifurcation scenario for zero through-flow - In the symmetry degenerate case without 
through-flow L-SPI (A ^ = 5, D > 0) and R-SPI (A = ^ 5, D < 0) have identical 
bifurcation properties. They are stable close to onset whereas RIB [A = B,D = {]) are 
initially unstable. In the driving range shown in flgures [T] and [2] the squared amplitudes of 
these two states grow practically linearly with the reduced distance /i from the common onset 
at /i = 0, albeit with different slopes. Then, there appear in a flnite supercritical driving 
interval stable CR-SPI solutions which transfer stability from SPI to RIB. The solution 
which bifurcates with B = Q out of the L-SPI is identifled as a L-CR-SPI with \A\ > \B\, 
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i.e., D > 0. The symmetry degenerate R-CR-SPI {\B\ > \A\,D < 0) bifurcates with A = 
out of the R-SPI. In the former \B\ grows and |y4| decreases - and vice versa in the latter 
- until the CR-SPI branches end with A = B,D = in the RIB solution. The amplitude 
variations of the CR-SPI solutions, however, are such that the sum S remains practically 
constant, cf. Fig. [2]^a). The RIB state loses stability outside the plot range of figures [H and 
[2] to another type of amplitude-modulated CR-SPI that are not discussed here. 

Through-flow induced changes - The axial through-flow significantly perturbs and changes 
structure, dynamics, and bifurcation behavior of the SPI vortex solutions discussed so far 



y, 0, 
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12l |: RIBs cease to exist in the strict sense, L- and R-SPI are no longer 



mirror images of each other, and also L-CR-SPI are no longer related to R-CR-SPI by this 
symmetry operation. However, spirals retain their spatiotemporal structure in the through- 
flow, i.e., they still do not depend on (p,z,t separately but only on the phase combination 
(pA = Mip + kz — ujAt or (pB = M(p — kz — cu^t, respectively. L-SPI (red color in Figs. [Tll2]) 
bifurcate for small Re > at a threshold value /ia < prior to R-SPI (orange color) which 
bifurcate at < /i^ ~ — /Ua out of the basic state. 

L-SPI are again stable at onset, but then lose stability to L-CR-SPI (violet color) which 
remain stable in the plotted parameter regime. The L-CR-SPI solution approaches the 
{Re = 0) RIB state with increasing /i, but retains with A ^ B a. finite distance D > 0. 
On the other hand, R-SPI are unstable for small and large /i, but stable for intermediate /x. 
Stability is exchanged with a R-CR-SPI (magenta color ) which has a stable as well as an 
unstable branch resulting from a saddle-node bifurcation at fis in Figs. [H O For small /i the 
unstable R-CR-SPI lies close to the [Re = 0) RIB solution and bifurcates with finite D < 
out of the R-SPI slightly above /i^. Note also that the sum of the squared amplitudes S of 
CR-SPI is no longer constant as for the case of Re = 0. 

Amplitude equations - The changes in spatiotemporal, bifurcation, and stability behavior 
of SPI, CR-SPI, and RIB states by a small through-flow can be explained and described 
close to onset quantitatively within an amplitude-equation approach. To demonstrate that 
we focus here on the bifurcation properties of the moduli \A\, \B\ of the critical modes. 

In order to reproduce the bifurcation and stability behavior of the aforementioned vortex 
states including the CR-SPI one needs coupled equations for A and B of at least quintic 
order. Higher-order terms that are suggested in l3| are not necessary to ensure the existence 
of CR-SPI solutions. Symmetry arguments [sl, m, Is] restrict the form of the equations for 
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the moduli to 

= \A\ [{fx - ixa) + bA\A\' + ca\B\' + CA {\A\'\B\' - \B\')] , (3a) 
= \B\ [(/i - /is) + bB\B\^ + cb\A\^ + CB (|i?| Vr - 1^1')] (3b) 

with real coefficients that depend in general on Re. Here we have discarded the quintic 
terms |y4|^ and \B\^ in view of the linear variation of the squared SPI moduli with /i, cf. 
Fig. [H Furthermore, we made a special choice for the coefficients of the last terms 
and I-BIIAI"^ in Eq. ([3]) that is motivated by the linear variation of the squared RIB moduli 
with /i for Re = and that suffices to describe the behavior for small Re as well. 

Coefficients - As a result of the inversion symmetry under z ^ —z which includes revert- 
ing the thro ugh- flow the coupled equations are invariant under the operation {A, B, Re) 
{B,A, —Re) so that the coefficients in Eq. ([3]) obey relations like, e.g., c^(i?e) = CB^—Re). 
To reproduce the bifurcation properties of the moduli for small Re as in Figs. [H [2] it suffices 
to incorporate the Re dependence to linear order in the coefficients 

fXAiRe) = -fx^^'^Re, fXB{Re) = fx^^^Re (4a) 
CA{Re) =c + c^^^Re, CB{Re) = c - c^^'^Re (4b) 

only and to ignore any Re dependence of the others by setting 

bA = bB = b, eA = cb = e. (4c) 

The choice bA = bs = b reflects the fact that the linear growth of the squared SPI moduli 
with the distance from their respective thresholds at fiA and is unchanged by the through- 
flow. On the other hand, the downwards shift of the L-SPI onset being for small through- 
flow of equal magnitude as the upwards shift of the R-SPI onset is reflected by HA^Re) = 
—fXBiRe) = —jx^^^Re with positive /i^^^. The coupling constant ca = = e ensures the 
existence of CR-SPI solutions [6]. The flow induced changes of the coupling constants ca 
and cb reflect the perturbation and destruction of the RIB states and ensures for positive 



Re L-CR-SPI solutions with D > when /i is large. The values of the coefficients [13|] were 
obtained from linear stability analyses and by fits to the full numerical nonlinear results. 
Fixed points - With the coefficients (jlj) it is straightforward to derive from ([3]) the following 
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relations for the fixed points of S and D 

IJ + {b + c)S- c^^^ReD - 2eD^ = 0, (5a) 
/i(^)i?e + {b- c)D + c^^^ReS + 265^ = 0. (5b) 

From these equations we have obtained the fines in tfie bifurcation plots of \A\'^ = S + D and 
= S — D in Fig. [Hand of S and D in Fig. [21 They compare very well with the symbols 
from the full numerical simulations of tfie Navier-Stokes equations, say, up to Re ~ 0.2 or so. 
Beyond tfiat fiigfier order terms in tfie i?e-expansion of tfie coefficients sfiould be included. 

Perturbation of the RIB CR-SPI bifurcation - Tfie black lines and symbols in Figs. [H 
[2] sfiow tfiat witfiout tfirougfi-flow tfie RIB state wfiicfi is stable at large ^ transfers stability 
to CR-SPl in a pitcfifork bifurcation located at ^2 in Fig. [H [2l Tfie symmetry and topology 
of tfiis unperturbed bifurcation out of tfie D = RIB solution at /i2 = — (& + c)S^:, S"* = 
(c— 6) / (2e) is best seen in Fig.[2](b). Tfiis figure sfiows also most clearly fiow tfie tfirougfi-flow 
witfi Re > perturbs tfiis bifurcation into a continuously varying stable solution brancfi 
witfi D > (violet) and a pair of stable and unstable solutions (magenta) tfiat are connected 
by a saddle-node bifurcation at fis- 

Tfie values of S and D at tfie unperturbed location /i = /i2 of tfie pitcfifork bifurcation 

5(/i2) = S2 = S, + s^Re^l^ + h.o.t., D{fi2) = = d2Re^''' + h.o.t., (6) 

grow oc Re^/^ and oc Re^^^, respectively, witfi constants = — (6+c)(/i*^^) + c^^^5'=i,)/(4e^) and 
S2 = 2ed\/ (6 + c). Tfie coordinates of tfie saddle-node bifurcation, as indicated in Fig. [2t^b) 
by tfie magenta arrows, 

/i5 = /i2 + f^sRe^^^ + h.o.t., Ds = dsRe^'^ + h.o.t., (7) 

vary oc Re^^^ and oc Re^^^ witfi tfie perturbation strengtfi Re. Here tfie constants are 
4 = [(c2 _ ^2)^(1) ^ 2e(6 + c)/i(i)] / (16e3) and fis = \{b + c)/i(i) - /iac^^) + Ae''d% \ /{2eds). 
Tfie comparison in Fig. [3] of tfiese lowest order approximations ([6], [7]) tfiat reflect an additive 
pitcfifork unfolding witfi tfie full fixed point result from Eq. sfiows good agreement for 
Re < 0.02. 

Larger through-flow - Here we investigate tfie transformation of tfie unperturbed RIB 
state into a CR-SPI as a function of tfie tfirougfi-flow strengtfi extending to larger Reynolds 
numbers up to Re = 4 wfiicfi are about 200 times larger tfian tfiose studied so far in tfiis 
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paper. For such large Re the amphtude equations with the simple Re dependence of the 
coefficients are no longer adequate to describe the behavior of the order parameters S 
and D. However, they continue to follow the same simple power laws with Re as the ones for 
small through-flow: When /i is well away from the location /i2 of the unperturbed pitchfork 
bifurcation as in Fig. Ka-c) one has S = + S^^'^Re^D = D^^'^Re + D^^^Re^. On the 
other hand, at /is, as in Fig. Kd-f) one has S = S^^^ + S^^'^\Re^'^f + 5(^/3)^^^1/3)4 ^ jj ^ 
/}(i/3)^gi/3 _|_ jjWi^jIq^/^Y. The lines in Fig. Hlshow fits of the data from the full numerical 
simulations (symbols) to these power law expansions for S and D. Fig. Ill^c,f) shows plots 
of the resulting squared moduli |Ap = S + D and \B\'^ = S — D. 

Note that inversion symmetry z ^ —z implies S{—Re) = S{Re) and D{—Re) = —D{Re). 
For Re > {Re < 0) the RIB state is transformed into a L-CR-SPI (R-CR-SPI) with A 
(B) being the major mode and B {A) the minor one. With increasing through-flow the 
major (minor) mode of the CR-SPI increases (decreases) until the minor mode vanishes in 
the transition to a pure SPI. 

Conclusion - A small applied flow is an inversion symmetry breaking perturbation for 
fluid wave structures that is easy to control. It significantly changes the bifurcation and 
stability transfer scenario involving symmetry degenerate TWs, mirror symmetric SWs, 
and mixed wave patterns. Numerical simulations of such waves in the counter-rotating 
Taylor- Couette system are explained and reproduced quantitatively by amplitude equations 
containing quintic order coupling and symmetry breaking perturbation terms. This approach 
can be applied to a wide range of 0(2) symmetric oscillatory instabilities in the presence of 
inversion symmetry breaking perturbations. 

This work was supported by the Deutsche Forschungsgemeinschaft. 
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FIG. 1: (Color) Bifurcation diagrams of L-SPI {Aj^O,B = 0,D> 0), R-SPI {A = 0,B ^0,D < 
0), L-CR-SPI {\A\ > \B\,D > 0), R-CR-SPI {\B\ > \A\,D < 0) and RIB (A=B, D=0 - only for 
Re = 0) as functions of fi and R2 for Ri = 240. Here, A and B are the amplitudes of the dominant 
modes U2,i(t) and U2^-i{t) ([1]), respectively, in the radial velocity field u at mid-gap. Solid (open) 
symbols denote stable (unstable) solutions of the full Navier-Stokes equations. Full (dashed) lines 
are the stable (unstable) solutions of the coupled amplitude equations. Red, orange, violet, and 
magenta show results for Re = 0.02. The color coding of the different solutions is given in the 
schematic 3D bifurcation diagram of (a). Therein the black dots denote Re = solutions without 
stability information. In (b) we show the projection of bifurcation diagrams onto the |yl|^ — 
plane. The color green refers to L-CR-SPI for Re = 0.1 where R-CR-SPI do not exist anymore. 
Black symbols and lines in (b)-(d) refer to Re = 0. 
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FIG. 2: (Color) Bifurcation behavior of the vortex states of Fig. [T] are shown here in plots of S (a) 
and D (b) versus fi and R2, respectively. The color coding is the same as in Fig. [TJ 
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FIG. 3: Power law behavior of quantities characterizing the break-up of the Re = pitchfork 
bifurcation from RIB to CR-SPI versus through-flow strength Re for Ri = 240. D2 and «S'2 refer to 
values at the location fi2 of the unperturbed bifurcation; Ds and ns are saddle node coordinates, 
cf. Fig. WIJo)- Solid lines come from the full fixed point equations ([5]); dashed lines refer to their 
leading order approximations (j6l7p . 
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FIG. 4: (Color online) Power-law variation of 5, D, |Ap = S + D, \B\^ = S — D with through-flow 
Re during the transformation from RIB (diamond) at Re = to CR-SPI (squares and triangles). 
The latter undergo a transition to SPI (circles) at large Re. Symbols denote results from the full 
Navier-Stokes equations for Ri = 240, lines show fits explained in the text. Right column: /x = /X2, 
i.e., the unperturbed pitch fork location; left column /j, = 0.011, i.e., far away from it. Note the 
different abscissa scaling of the two columns. 



